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We calculate the disconnected contribution to the form factor for the semileptonic decay of a 
D-meson into a final state, containing a flavor singlet eta meson. We use QCDSF «^ = 2 + 1 
configurations at the flavor symmetric point /«„ — m^ = m^ and the partially quenched approxi- 
mation for the relativistic charm quark. Several acceleration and noise reduction techniques for 
the stochastic estimation of the disconnected loop are tested. 
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1. Introduction 



Semileptonic decays of D-mesons contain rich physics. Lattice calculations of the form factors 
for these decays are important for the search for hints of new physics through the determination of 
CKM matrix elements. These form factors have been well-studied on the lattice. Previously, some 
of us tested a stochastic method to measure 3-point functions needed to calculate the semileptonic 
decay form factor [|l|]. The advantage of stochastic methods is that we have access to a greater range 
of momenta at fixed cost. This enables us to extract the form factor more reliably from results for 
the three point functions at different momentum transfers. 

In particular, the Dg meson is interesting for flavor physics. Its major semi-leptonic decay is to 
77 and Tj', which has a contribution from a disconnected loop diagram (Fig. [l|). The loop runs over 
three light flavors so the effect is enhanced by a factor three, and thus may be large. The purpose 
of this work is to test the feasibility of measuring the disconnected diagram, and to quantify its 
contribution to the form factor. 

We extract the scalar form factor /o from the relation [^] : 

Mq^) = ^ ji^lSlDs), (1.1) 

mj)^ — 

where S = /c is a scalar current made from charm and light quarks, m, are the masses of the quarks 
and mesons. The matrix element can be extracted from the following ratio of 3-point over 2-point 
functions: 

{ri(k,ti)\S{q,t)\D,{p,tf)) =Zr,ZD, ^'^^^~^'f ~^f'^''^\ =Z^ZDM.q,P.t-ti,tf-ti), (1.2) 

C{{t -ti;k)C2"{tf -t;p) 

and similarly for r]' . For large tf — ti and t — ?, this ratio should approach a constant. Z-q and Zu, are 
the overlap factors between the meson state and the interpolating operator, which can be extracted 
from the two point functions C2 {t — tfX) and C^'(f/ — t; p), respectively. The two point functions 
for T] and r]' also have a disconnected part, however, at mps ~ 445 MeV we expect its contribution 
to the mass to be small and we neglect it in this first exploratory study. 

We use QCDSF 24^ x 48 «/ = 2 + 1 configurations [|3|]. So far we only use the SU(3) symmet- 
ric set {Ki = Ks = 0.1209) with lattice spacing a ~ 0.08 fm. This was generated using the tree-level 
Symanzik-improved gluonic action and non-perturbatively improved Wilson fermions with stout 
links in the derivative terms (SLiNC action). We use the same relativistic quark action for the 
(quenched) charm quark with Kcbavm = 0.11. Note that since we use the flavor SU(3) symmetric 
configurations, the disconnected contributions in the — )• T] 3-point function cancel, when we 
identify r\ = r\^. The Chroma software package [Q] is used for some of the analysis. 
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Figure 1: Connected (left) and disconnected (right) diagrams which contribute to Ci{tf,t;p,q). We use a 
stochastic method to estimate the all-to-all propagators, denoted by blue lines. 
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2. Noise Reduction techniques 

In order to calculate the disconnected loop, all-to-all propagators are required. These are 
estimated using stochastic methods, which involve performing A'^ inversions of the light quark Dirac 
operator for each configuration; should be large enough to give sufficiently small stochastic 
errors relative to the gauge noise. For some quantities the stochastic noise dominates the overall 
uncertainty and it is important to use efficient noise reduction techniques. 

We measure the disconnected "loop" 

where M is the Dirac operator for a light quark and is a smearing function. The stochastic 
estimation of the all-to-all propagator M^'(xV;-^",0 involves the following approximation: 

A/-'4fwtol + ^(^). (2.1) 

where It/,) is a random noise vector and \si) = M^^Itj,). We use -^(^2 + '^2) complex random 
numbers for the noise vector. For each / we need to smear both |t],) and (|t],) must be smeared 
after solving for \si)) so we need 2N applications of the smearing operator. This significantly 
increases the computer time needed to calculate the disconnected loop. Time dilution (partitioning) 
is implemented: the noise vector is only non-zero on one or two time slices. 
We test the following three noise reduction techniques. 

Spin dilution/partitioning This uses projected noise vectors on a single spinor component and 
sums over the projections afterwards [^: 

^it\^l"Mvt\ (2-2) 

■'^ a=li=l 

where {rj-"'^) = P^^^rji) is the projected noise vector. It requires 4N inversions but for some 
quantities the stochastic error is reduced by a factor greater than 2. In addition, we can 
reduce the cost of smearing because the spin projection P^"^ commutes with the smearing of 
our choice. A naive scaling gives SA'^ smearing operations, but we only need 5A^ applications: 
4A'' for Is^"^) and A'' for |t],-). 

Hopping Parameter Acceleration (HPA) [^] This is based on the following identity 

(k-D)"M"' = M-'^ - kD - (kD)^ {kD)"~\ (2.3) 

where kD is the hopping part of the Dirac operator. Note that the derivative operator satisfies 
tr[75K"D] = due to the spinor structure so that this term only contributes to the noise. This 
means that (k"D)^M^' represents an improved estimate of M^^ (we call it n = 2 HPA). As 
long as the smearing is diagonal in spinor space, this is also true for the smeared all-to-all 
propagator. 
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Truncated Solver Method (TSM) For some quantities the ultra violet modes dominate. In these 
cases, using a small number of CG iterations in the solver for the solution vector \si) provides 
a good approximation, for example, to the disconnected loop [Q, ^J. To arrive at an unbiased 
estimate, a correction term needs to be added to the truncated part: 

Y Ni y N1+N2 

= TrY.\^tmnc,i){rii\ + — £ kbias,;>(r7;l- (2-4) 

■'^1 1=1 ■'^2 j=Ni+l 

The first term uses the truncated solution |i'trunc,i)> which is cheap to calculate and typically 
causes the main part of the stochastic error. The second term contains l^bias,^) = kconv,;) — 
ktrunc,;)> whcrc |i'conv,;) is a Converged solution. I^'conv,;) is expensive, and only accounts for a 
small part of the stochastic error if Isbias,;) does not contribute significantly to the observable. 
Therefore, by tuning parameters — n: number of CG-iterations for the truncated part, A'^i: 
number of stochastic noises for the truncated part, A'2: number of stochastic noises for the 
bias part — we can reduce the total calculation cost. We use a CG solver for the truncated 
solutions and a BiCGstab solver for the converged solutions. 



3. Comparisons 

We investigate the noise reduction techniques using one configuration. We use Wuppertal 
smearing [^] for the quarks, with parameters which are tuned to minimize the contributions from 
the excited states to the effective mass. 

In Figs. ^ and ^ we plot the stochastic errors for various combinations of the noise reduction 
techniques. In each case, the computational cost is fixed. The horizontal axes correspond to n, the 
number of iterations of the solver in the TSM. The data at « = — 100 indicate the results without 
the TSM. In particular, the red plus symbols ("+") show the results without any noise reduction 
techniques. For a fixed n, we have optimized A^i and to give the smallest stochastic error under 
the cost condition 

Nl (WTCG + Tsmear) + N2 (ntcG + «conv ^BiCGstab + ^smeai-) = Constant, (3.1) 

assuming the square of error, o^^^^^^ , to scale according to 



c^s1och. = tf + tf, (3-2) 



-2 — 

Ni N2' 

where /i and fx are the variances of the first and second terms in eq. ( |2.4[ ), respectively. ?iconv is the 
number of iterations needed to obtain the converged solution. Tcgi 'J'BiCGStab and Tgmear represent 
the computer time needed for 1 CG iteration, 1 BiCGstab iteration, and smearing, respectively. The 
optimal ratios of N1/N2 are around 1 (10), with (without) smearing. 

Although small differences between the results are not significant due to the uncertainty on the 
stochastic errors, in all cases spin dilution together with HPA (purple squares), gives the minimum 
error when combined with TSM. Therefore we use this combination in the following analysis. 

The gain factor, 

(without noise reduction) 



a^(with noise reduction) 



(3.3) 
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TSM: optimal stochastic errors, with fixed cost : pOOO, no smearing 



TSM: optimal stochastic errors, with fixed cost : pOOO, with smearing 
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Figure 2: Estimated stochastic errors at fixed cost for p 
Data at n = —100 are without TSM. Left panel: without 
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Figure 3: The same as Fig. but for p = (1,0,0). 



strongly depends on the smearing. Without smearing (left panels), we obtain maximum gain factors 
of 16 - 25, which translates into a reduction of the computational cost of the same magnitude. With 
smearing, it is only about a factor 2. This is because the contribution to the error from the bias part 
(i.e., /2) is larger than or of the same magnitude as the truncated part (/i). 

4. Results 

Having optimized the noise reduction, we can now measure the disconnected contribution to 
the form factor. For the TSM, we truncate after n = 20 CG iterations and the numbers of noise 
vectors are Ni = 10 and = 20. A total of 939 configurations were used in the analysis. 

Following our previous study [Q], we use stochastic techniques for the connected contribution 
as well. The noise vectors are placed at the sink of the D., meson (denoted by a red circle in Fig. [I]). 
For each configuration, 24 x 4 spin diluted noise vectors were computed for the charm quark. In 
terms of momenta, 57 different combinations of p for the Dg meson were calculated. Note that a 
similar calculation with the sequential method would require 57 x 12 inversions. 

In order to extract the matrix elements in eq. ( |1.2| ), we fixed the time separation between the rj 
source and the sink separately for the connected (?/ = 24, ti = 0) and the disconnected (?/ = 24, 
ti = 16) matrix elements. We combine the two contributions afterwards. For the connected part, 
taking the maximum separation tf — ti = T /2 = 24 enables us to average over the forward and 
backward propagations. For the disconnected part, in order to average the forward and backward 
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Figure 4: Ratios of 3-point over 2-point functions, R in eq. 1.2, for connected and disconnected parts. 
k^q~ (0,0,0) for the left panel and ~k — q~ (1,0,0) for the right panel. 



propagations, the noise vector has a non-zero value at two time-slices separated by 16 time-slices 
(f/±8). The usage of different tf — tj for the connected and disconnected 3-point functions is 
allowed because we have assumed = m,,/ (remember that »!„ = = nis). 

Fig. ^ shows the ratio of the correlation functions, which corresponds to /^{q^) /Zr^Zr,^. The 
disconnected part is multiplied by 3 because of the 3 light flavors. The errors for the disconnected 
contribution are small enough to obtain signals, significantly different from zero. 

In Fig. |5| we show the form factors for the octet {r]^) and singlet (tji) T]s: 

Ir/g) = -^{\uu) + \dd) — 2\ss)) connected only, (4.1) 

v6 

Itji) = -^{\uu) + \dd) + \ss)) connected — 3 X disconnected. (4.2) 

v3 

Preliminary fits to foiq^) of the form f(){q^) = give /o(0) = 0.75(3) and /o(0) = 0.52(5), for 
Ds — )• T]8 and — )• t]i, respectively. Also included in Fig. |5| is a value from light cone QCD sum 
rules for the decay into Tj [|lO|], /o(0) = 0.45(14). Due to SU(3) flavor symmetry the — ^ r\?, form 
factor also represents the form factor of D — Ivn and D — IvK. Note that /o(0) for r]\ is smaller 
than that for r]^. This is consistent with the form factors for B — tj, T]' [|l[], which is the heavy 
quark limit, can be compared to our calculation. 



5. Conclusions 

We tested three methods (and their combinations) of noise reduction techniques for measuring 
the disconnected contributions to the meson semi-leptonic decay form factor. The combination 
of spin dilution, hopping parameter acceleration and truncated solver method was found to give the 
biggest gain in computer time. These noise reduction techniques allowed us to measure non-zero 
contributions to the form factor, on SU(3) flavor symmetric QCDSF Uf = 2 + \ configurations. 
Further studies with non-SU(3) symmetric Uf = 2 + \ configurations are planned. 

This work was supported by the EU ITN STRONGnet (grant number 238353) and the DFG 
SFB/Transregio 55. SC acknowledges support from the Claussen-Simon-Foundation (Stifterver- 
band fiir die Deutsche Wissenschaft). JZ is supported by the Australian Research Council under 
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Figure 5: Form factor foiq^) for — > ZvTjg and — > Ivrji. Errors are statistical only. A value from QCD 
light cone sum rules (LCSR) | |l0[ | is also plotted. To enhance visibility, the fitted values at <7^ = and the 
LCSR results are slightly shifted to the right and left, respectively. 

grant FT 100 100005. The calculations were performed on the Athene HPC cluster at the University 
of Regensburg. 

References 

[1] R. Evans, G. Bali and S. Collins, Improved semileptonic form factor calculations in lattice QCD, 
Phys. Rev. D 82 (2010) 094501 [arXiv: 1008.3293 [hep-lat]]. 

[2] H. Na, C. T. H. Davies, E. Follana, P. Lepage and J. Shigemitsu, D semi-leptonic decay form factors 
with HISQ charm and light quarks, PoS LAT2009 (2009) 247 [arXiv:0910.3919 [hep-lat]]. 

[3] W. Bietenholz et.al.. Flavour blindness and patterns of flavour symmetry breaking in lattice 
simulations of up, down and strange quarks, Phys. Rev. D84 (201 1) 054509 [arXiv:l 102.5300 
[hep-lat]]. 

[4] R. G. Edwards and B. Joo [SciDAC, LHPC and UKQCD Collaborations], The Chroma software 
system for lattice QCD, Nucl. Phys. Proc. Suppl. 140 (2005) 832 [arXiv:hep-lat/0409003]. 

[5] S. Bernardson, P. McCarty and C. Thron, Monte Carlo methods for estimating linear combinations of 
inverse matrix entries in lattice QCD, Comput. Phys. Commun. 78 (1993) 256. 

[6] C. Thron, S. J. Dong, K. F. Liu and H. P. Ying, Pade-Z(2) estimator of determinants, Phys. Rev. D 57 
(1998) 1642 [arXiv:hep-lat/9707001]. 

[7] S. Collins, G. Bali and A. Schafer, Disconnected contributions to hadronic structure: a new method 
for stochastic noise reduction, PoS LAT2007 (2007) 141 [arXiv:0709.3217 [hep-lat]]. 

[8] G. S. Bali, S. Collins and A. Schafer, Effective noise reduction techniques for disconnected loops in 
Lattice QCD, Comput. Phys. Commun. 181 (2010) 1570 [arXiv:0910.3970 [hep-lat]]. 

[9] S. Giisken, U. Low, K. H. Mutter, R. Sommer, A. Patel and K. Schilling, Nonsinglet axial vector 
couplings of the baryon octet in lattice QCD, Phys. Lett. B227 (1989) 266. 

[10] K. Azizi, R. Khosravi and F. Falahati, Exclusive D^ — ?> (tj , ri')lv decays in light cone QCD, J. Phys. G 
38 (2011)095001 [arXiv: 101 1.6046 [hep-ph]]. 

[11] R Ball, G. W. Jones, B r\'''\form factors in QCD, JHEP 0708 (2007) 025 [arXiv:0706.3628 
[hep-ph]]. 



I. 



fit ■ 

Ti, LCSR I — 



7 



